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ABSTRACT 

We propose that the very narrow (FWHM<10keV) gamma-ray lines observed at 
1.495 and 6.695 MeV from SS433 (Lamb et al. 1983, 1984) are blueshifted 1.369 
and 6.129 emissions from deexcitations of 24 Mg* and 16 0* in grains moving with 
tne jets and i nel asti cal ly excited by interactions with the ambient medium. 

As previously shown (Lingenfelter and Ramaty 1977), energetic particle 
interactions in grains produce very narrow gamma-ray lines from deexcitation 
of nuclear levels whose lifetimes are long enough that the excited nuclei stop 
before deexcitation. The presence of grains in the jets resolves hitherto 
discussed difficulties of inelastic-excitation models for gamma-ray production 
in SS433, namely the very narrow widths of the observed lines and the absence 
of otner strong lines expected from abundant elements. The model that we 
propose could be distinguished from a previously proposed fusion model (Boyd 
et al . 1984) Dy gamma-ray line observations. 


*In press, Astrophys. J. Letters, 1984. 
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I. introduction 

Lamb et al . (1983, 1984) have observed blue- and red-shifted gamma-ray 
lines from SS433 with the HEAO-3 high-resolution gamma-ray spectrometer. The 
lines with the highest statistical significance, observed at 1.495 and 6.695 
MeV, are very narrow (FWHM < 10 keV) and their central energies are determined 
to an accuracy of a few keV (Lamb et al. 1983, 1984). Boyd et al. (1984) have 
interpreted these lines as blueshifted 15 0* deexcitation lines of rest 
eneryies 1.380 MeV and 6.175 MeV. These two lines result from the cascade 
deexcitation of the 7.556 MeV level of 15 0, through its 6.176 MeV level, to 
tne ground state. The 15 0 level at 7.556 can be populated by thermonuclear 
fusion of l4 N and protons through a narrow resonance at 0.278 MeV. The ratios 
of tne observed line energies, 1.495/6.695 = 0.2233, and line intensities, 
<f(6.695)/<p( 1 .495) '■0.7, are, within observational uncertainty, the same as 
those expected from the l5 0 lines, 0.2235 and 1, respectively (Ajzenberg- 
Selove 1981). 

There is, however, an alternative identification of the observed lines. 
Lamb et al. (1983) orginally identified the 1.495 MeV line as a blueshifted 
1.369 MeV line resulting from 24 Mg*, excited by inelastic collisions with 
protons. One of the strongest lines that is expected to accompany this line 
is that at 6.129 MeV from 16 0* (Ramaty, Kozlovsky and Lingenfelter 1979). By 
a remarkable coincidence, the ratio of these two line energies, 0.2233, is 
almost exactly the same as that of the two 15 0 lines and is, within 
observational errors, precisely equal to the ratio of the energies of the two 
observed blueshifted lines. Thus, we propose that the observed lines from 
SS433 are due to 24 Mg* and 16 0* in the jets, excited by inelastic interactions 
with ambient protons. 

In both the inelastic-excitation and the fusion models, the implied 
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Dluesnift is consistent with that expected at the epoch of observation from 
the approaching jet of SS433. As suggested by Lamb et al. (1983), a weaker 
line, observed at approximately 1.2 MeV, could then be the redshifted 
counterpart of the 1.495 MeV line produced in the receding jet. For a recent 
discussion of the kinematics of SS433, see Margon (1984). 

The lines observed from SS433, however, cannot be produced by excitations 
of 24 Mg and 16 0 in gas. For at proton energies above the line production 
threshold, recoils of the l6 0* nuclei broaden the 6.129 MeV line to a width of 
- 100 keV which is totally inconsistent with the observed width (< 10 keV). 
Moreover, at such bombarding energies and for galactic abundances (Cameron 
1982), the two strongest lines would be at 6.129 MeV and 4.438 MeV, and not at 
1.369 MeV as observed from SS433. 

Both of these difficulties, however, can be resolved if the gamma-ray 
lines are produced in grains or other dense particles (p - 1 g cm" 3 ) moving 
witn the jet velocity. Lingenfelter and Ramaty (1977) showed that energetic 
particle interactions in grains can produce very narrow gamma-ray lines, and 
that strong grain lines are expected at 6.129 and 1.369 MeV from 16 0 and 24 Mg 
deexcitation. These gamma-ray lines are very narrow because the lifetimes of 
the corresponding excited levels are sufficiently long (> 1 psec) that the 
recoiling nuclei stop before they emit gamma rays. On the other hand, the 
lifetime of the 4.439 MeV level of 12 C is too short to allow the excited 12 C 
nuclei to stop. Therefore, the recoi 1-broadened 4.438 MeV line would be so 
wide (FWHM ~ 100 keV) that it would fall below the instrumental sensitivity of 
tne HEAO-3 detector. Furthermore, the well-known depletion of Me and 0 in 
grains with respect to Mg, Si and Fe could explain the observed relative line 
intensities without having to assume arbitrary elemental abundances. Overall 
tnen, the assumption of the existence of grains in the jets allows us to use 



the simplest interaction model, namely that the nuclear excitations result 
from the interaction of the jets with the ambient medium. 

II. VERY NARROW LINES FROM GRAINS 

We now consider the strongest very narrow lines expected from excited 
nuclei that lose tneir recoil energy in dense particles before deexcitation. 
The size of these particles, which we will hereafter refer to simply as 
grains, must be less than ~ 1 cm to allow them to be transparent to gamma rays 
and greater than ~ 10 -4 cm to stop the excited nuclei. The principal grain 
lines are listed in Table 1, together with the gamma-ray emission processes, 
mean lifetimes and principal excitation processes (Ramaty, Kozlovsky and 
Linyenfelter 1979). 

The relative intensities of the very narrow lines depend on the 
excitation cross sections, the elemental abundances and the probability that 
the recoiling nucleus will stop before emitting gamma rays. The cross 
sections, a, given in Table 1 are for an interaction energy of 33 MeV, 
corresponding to tne velocity (V=0.260c, Margon 1984) of the ambient protons 
in the rest frame of the jets. These cross sections are from Ramaty, 

Kozlovsky and lingenfelter (1979), with corrections inferred from the 
measurements of Dyer et al. (1981) at energies less than 23 MeV. 

We assume that essentially all the heavier, less volatile elements, Mg, 
Si, Fe, etc., are in the grains with the relative abundances of Cameron 
(1982), and that the more volatile lighter elements are depleted, such that in 
the grains C, H and 0 are reduced by a factor f relative to the Cameron (1982) 
abundances and there is essentially no H, He, Ne or S. Such abundances (Table 
1) are consistent with studies (e.g. Morton 1975) of gas depletion in the 
interstellar medium. We treat the factor f as a free parameter, which we 
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expect to be about 1/2 or less. 

The probabi 1 i ties that the recoiling nuclei stop in the grains before 
emitting gamma rays were calculated in detail (Ramaty, Kozlovsky and 
Lingenfelter 1979). For grains larger than a few times 10" 4 cm, these 
probabilities are essentially unity for all the lines listed in Table 1, 
except the 4.438 MeV line for which it is vanishingly small. 

From these cross sections, abundances and stopping probabilities, we 
calculate the relative intensities of the grain lines listed in Table 1. As 
can be seen, the two strongest lines are at 6.129 and 1.369 MeV. The observed 
ratio, q>( 6.695) /<p( 1 .495) - 0 . 7 , implies that the 16 0/ 24 Mg ratio in the grains 
is ~ 4, or that f - 0.17. Such a fraction is consistent with the locking of 0 
in grains due to the formation of chemical compounds with Mg, Si and Fe. 

In addition to the observed lines, the inelastic-excitation model also 
predicts lines at 0.847, 0.931, 1.317, 1.634 and 1.779 MeV with intensities 
about 1/2 tnat of the 1.369 MeV line. The relative intensities of the 0.847 
and 1.779 MeV lines, however, depend on the Fe and Si abundances. But, 
independent of the abundances, there should be 1.634 MeV emission from Ne 
produced by 24 Mg spallation. Norman and Bodansky (1984) have suggested that 
the intensity of the 1.634 line relative to that of the 1.369 MeV line should 
be ''0.85, based on the cross sections of the reactions 24 Mg(p ,p‘ ) 24 Mg* and 
24 My (p ,x) 20 Ne* at 33 MeV, and that this high ratio would seem to rule out 
inelastic-excitation at such an energy since the 1.634 MeV line was apparently 
not seen. However, as we show in Table 1, the ratio of these lines depends 
also on the abundance of 28 Si , because of the strong contribution of 28 Si 
spallation to the 1.369 MeV line. Therefore, the apparent absence of the 
1.634 MeV line is not an argument against the inelastic-excitation model. 

Nuclear interactions in 0, Mg, Si and Fe also produce positrons. We 
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estimate that at 33 MeV about 1.8 positrons should be produced for each 1.369 
MeV photon. The width of the 0.511 MeV line, however, will depend on the 
temperature and density of the annihilation site. 

III. ABSOLUTE LINE LUMINOSITIES, ENERGY LOSS, GRAIN TEMPERATURE AND X-RAY 
PRODUCTION 

We now investiyate the energetics of the gamma-ray line emission. 
Considering for simplicity the instantaneous interaction rate between the jet 
grains and the ambient gas, the expected luminosity in the very narrow 1.369 
MeV line is given by 

L( 1.369 MeV) = »„ g * N $1 ^ ] V E^ (1) 

where N^g and N 5 ^ are the total numbers of Mg and Si nuclei in the grains, n H 
is ambient hydrogen density, and E y = 1.369 MeV. Using the values of a and 
the abundances given in Table 1, we obtain 

L( 1 .369 MeV) * 3.3xl0 34 n H ( M gr / M Q ) erg/sec, (2) 

where Mgp/My is the total mass of the grains in solar masses and n H is in 
cm -3 . This nuclear deexcitation radiation is accompanied by Coulomb energy 
loss, bremsstrahlung and a variety of other radiations which provide further 
constraints on the model. 

The Coulomb energy loss, W c , is due to ambient protons and electrons 
traversing the grains and gas in the jets with energies of 33 MeV and 18 keV, 
respectively, in the jet rest frame. Taking the bulk of the jet gas to be 
hydrogen of mass M gas , we obtain 
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where the four expression for dE/dx are energy loss rates per g cm -2 of 
protons and electrons in grain material and hydrogen. From Barkas and Berger 
(1964) and Berger and Seltzer (1964) we find 


W c - 6.7xl0 38 n H (M gr + 2.6 M ga$ )/M 0 erg/sec. (4) 

Combining equations (2) and (4), and substituting the observed average 

bluesnifted 1.369 MeV luminosity (~ 10 37 erg/sec, Lamb et al. 1983), we obtain 

fl 3 2xl0 41 (1 + 2.6 M /M ) erg/sec. This large Coulomb energy loss must be 
c s 9 r 

compensated by continuous acceleration, since otherwise the jets would slow 
down and the Doppler shift would change with time, contrary to the observation 
of an essentially constant blueshift and very narrow line for a period of 
several days (Lamb et al. 1984). Such constancy of the jet velocity may be 
maintained in condensed material by line-locking (Pekarevich, Piran and Shaham 
1984). 

The minimum = 2 x 10 41 erg/sec, obtained when M gas << M gr , exceeds the 

Eddington luminosity of a solar mass compact object by about 3 orders of 

magnitude. Therefore, should be released primarily as mass motion in the 

jets. Letting the Coulomb energy loss equal the jet kinetic energy flow, 

*c * oy + tigas^ y2 /2, and assuming that l*l gr /fl gas = M gr /M gas , we find 

that M 3 10" 4 e M n /y , where e is a numerical factor between 1 and 2.6, 
gr 0 

depending on M gas /M gr . 

We estimate that the Coulomb energy loss in the grains will heat them to 
a temperature high enough to drive off the more volatile compounds but not the 
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most refractory materials, particularly MgO. The equilibrium temperature of 
the grains, T, can be estimated by letting the rate at which energy is 
deposited in them, 2 x 10 41 erg/sec, equal a T 4 A, where A is the total surface 
area of the grains. This area can be calculated from the mass of grains 
involved in nuclear excitation at any one time and the area-to-mass ratio of 
the grains, A/M = 3/pr - 10 4 cm 2 /gm, assuming a density p ~ 3 gm/cm 2 and a 
radius r - 10“ 4 cm. Such a size and density are large enough that the thermal 
emission will be essentially blackbody and the bulk of the recoiling, excited 
My and 0 will stop in the grains to emit very narrow gamma-ray lines. The 
mass of grains in the emission region depends on ftg r > estimated above, and the 
length of the emission column, which could be 4 days travel distance or more, 
based on the length of flaring (Seaquist et al . 1982) If we take - 4 days 
and M yr - 2.6 x 10’ 4 M 0 /yr, we find a mass of ~ 3 x 10" 6 M 0 corresponding to a 
radiating area of 6 x 10 31 cm 2 and an equilibrium temperature of - 2750°K, 
significantly below the 3125°K melting temperature of MgO. We emphasize, 
however, that the relative line intensities given in Table 1 are for grains 
that consist not only of Mg and 0 but also C, Si and Fe. 

The ambient electrons and protons traversing the grains and gas in the 
jets produce a bremsstrahlung luminosity < 18 keV of 


L.. = 2 Vn 


[ (££) M + 
L 'cTjT'r,gr gr 


(dE\ 
Wr, 


H M gas ) ’Z*U> 35 " f 


("gr * °' 3 "gas> /M 8 


erg/sec.(5) 


Comoining equations (2) and (5) and substituting the observed line luminosity, - 

we find L (< 18 keV) - lxlO 38 (1 + 0.3 M /M ) erg/sec. These X rays must 
x ~ gas gr 3 J 

be substantially attenuated since the minimum bremsstrahlung of 10 38 erg/sec 
exceeds the observed < 18 keV luminosity (~ 10 36 erg/sec, Marshall et al. 

1979) by about a factor of 100. 
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The X rays, as well as the the gamma rays, could be attenuated or 
aosoroed in both tne ambient medium and the jets. We first estimate the 
attenuation in the ambient medium. We assume that the characteri sti c path 
length is comparable to tne length of the gamma-ray emitting portion of the 
jets, which we approximate by i = V Mgp/'^gp* From equation (2), with 
L( 1.369 MeV) = 10 37 erg/sec and = ICf 4 Mg/year, we find a column 
deptn i n u = 7xl0 23 cm" 2 . For such a value, the ambient medium is quite 

n 

transparent to 1.5 MeV gamma rays, x - 0.1, but for X rays, it has an optical 
depth of x = l at ^ 10 keV for Cameron (1982) abundances, and it could be much 
more opaque if the heavy element abundances were enhanced. 

We now estimate the opacity of the jets assuming cylindrical geometry 
sucn tnat t = anR = atf/xAVR, where R is the radius of the jet, A is the mean 
atomic mass and a is the mean absorption cross section (Plechaty, Cullen and 

Howerton 1975). For the minimum opacity, correspondi ng to M gas << Mg r 
and M - 10" 4 Mg/y, w e obtain optical depths across the jets 
of x = 1.4xlO lo /R at 1.5 Me'/ and x = 1.8xl0 13 /R at 10 keV, where R is the 
radius in cm. Thus, for the jets to be transparent to gamma rays, R must 
exceed ~ 10 LO cm. At such a value of R, the X-ray opacity is so large that 
tne Dremsstrahlung would only be ~ 10% of that observed. Thus, the expected 
X-ray oremsstrahlung is consistent with observations. 

IV. CONCLUSIONS 

We have studied the o-rigin of the recent HEA0-3 gamma-ray line 
ooservations of SS433. Although no independent confirmation of these 
ooservations has yet oeen reported, the correlation of the gamma-ray emission 
with radio flaring and the agreement of the gamma-ray source direction with 
tne known direction of SS433 are good arguments that the observed lines (Lamb 
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et a 1 . 1983, 1984) are indeed from SS433. Strong additional support for this 
identification is provided Dy the fact that a single jet velocity, appropriate 
to tne epocn of the gamma-ray observation, can accurately account to a 
precision of a few tenths of a precent for the line energies of the two 
observed bluesnifted lines in terms of rest energies predicted by theory. 

We suggest that tne gamma-ray line emission results from the inelastic 
excitation of 24 Mg and 16 0 in grains moving witn the velocity of the jets and 
interacting with the ambient gas in the acceleration region. This model 
provides an alternative identification of the observed lines to that of the p- 
14 N fusion proposed by Boyd et al. (1984). Since both models predict 
additional gamma-ray lines at differing energies, observations should clearly 
distinguish between them. 

The essential ingredient of the inelastic-excitation model is the 
presence of dense particles in the jets. The possibility of highly clumped or 
condensed matter in tne jets was also suggested recently by Pekarevich, Pi ran 
and Shaham (1984), if line locking is responsible for the accel eration. Much 
further study is needed, however, to show how these particles are produced, 
now long they survive in the jet environment, and what role they might have in 
the unsolved acceleration mechanism. 
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TABLE 1 






VERY 

NARROW LINE EMISSION 



Pnoton 

Emi ssion 

Mean 

Exci tati on 

o(mb) 

Abundance 

Li ne FI ux 

Energy 

Process 

Li fe 

Process 



(MeV) 


(psec) 





0.847 

56p e *0.3C7^ g>s 

9.7 

56 Fe(p,p‘ ) 56 Fe* 

150 

1.0 

0.5 

0.931 

56 Fe *0.93 W g>s 

13 

ss Fe(p,pn) 55 Fe* 

180 

1.0 

0.6 

1.317 

56p e *1.317^g #s 

' 1.4 

56 Fe(p,pn) 55 Fe* 

150 

1.0 

0.5 

1.369 

24 My* 1 • 3 69 -*-g.s . 

1.8 

^Mg(p ,p ‘ ) 24 Mg* 

130 

1.0 





23 Si (p ,x ) 2I+ Mg* 

90 

1.1 

1.0 

1.634 

20 Ne*I* 63 ‘ + -^. S . 

1.2 

20 Ne(p,p‘ ) 20 Ne* 
24 Mg(p,x) 20 Ne* 

90 

150 

- 0 
1 

0.5 

1 .779 

2S$i * 1. 7? 9-^g >s t 

0.7 

28 Si ( p , p 1 ) 23 Si* 

90 

1.1 

0.4 

4.438 



>2C(p,p')‘ 2 C* 

16 0(p,x) 12 C» 

70 

13f 




C .06 

80 

2 2 f 

- 0 

6.129 

U 0 *6.13U g>s< 

24 

1 6 0( p ,p 1 ) 16 0* 

20 Ne(p,x) 1S 0* 

50 

22f 




60 

- 0 
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